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NATIONAL ADVISOR! COMMITTEE FOR AERONAUTICS 


TECHNICAL MEMORANDUM 1234 


COMPUTATION OF THIN-WAILED PRISMATIC SHELLS* 
By V. Z. Vlasov 


1. Fundamental Assumptions 

We consider a prismatic shell consisting of a finite number of 
narrow rectangular plates and having in the cross-section a finite 
number of closed contours (fig. 1(a)). We shall assume that the 
rectangular plates composing the shell are rigidly joined so that there 
is no motion of any kind of cne plate relative to the others meeting 
at a given connecting line. The position of a point on the middle 
prismatic surface is considered to be defined by the coordinate z, 
the distance to a certain initial cross-section z = 0, and the 
coordinate s determining its position on the contour of the cross- 
section. 

Let the function u(z,s) represent a longitudinal displacement 
of the point (z,s), that is, a displacement in the direction of the 
generator (positive in the direction of increasing z) and the 
function v(z,s) represent a tangential displacement, that is, a 
displacement in the direction of the tangent to the contour of the 
cross-section (positive in the direction of increasing s). 

We shall represent these displacements in the form of the sums: 


u(z,s) 


m 

U i (z)(p i (s) 


(i = 1,2,..., m) (1.1) 


v(z,s) 


n 

> V k (z>|r k (s) 


(k = 1,2,..., n) (1.2) 


■where the functions U^(z) and V k (z) depending only on z are the 
required functions and the functions <pp(s) and \[r k (s) are subject to 

a preliminary choice. The magnitudes m and n will be explained 
below. 
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We consider an elementary tranverse lamina cut out "by the planeB 
z = constant and z + dz = constant. Such a lamina will he regarded as 
a plane system of rods, that is, a frame consisting of closed contours. 

Let us consider the deformed state of this lamina determined only 
by the longitudinal displacements u(z,s). The plane contour of such an 
elementary frame, while remaining on the prismatic surface of the shell, 
goes over into a space line that can be determined relative to the 
initial section z = constant for the chosen functions cp^(e) by equa- 
tion (l.l). We make the assumption that, in the case of the deformed 
state determined only by the longitudinal displacements, the rectilinear 
elements of the frame, starting from the plane z = constant, remain 
straight. Such an assumption is equivalent to the hypothesis of plane 
sections assumed individually for each of the narrow rectangular plates 
composing the given shell. In this case the position of an elementary 
frame after deformation is entirely determined by the longitudinal 
displacements of its m Joints relative to the plane z = constant. 

Hence the elementary frame can be considered as a rod system possessing 
relative to the longitudinal displacements m degrees of freedom. 

We shall take in equation (l.l) the required functions U-i(z), 

Ug(z), ..., ^n(z) as the longitudinal displacements of the m Joints 
of an elementary frame. The functions q>p(e), <P 2 ( s )j •••> ^(s) 

corresponding according to equation (l.l) to these displacements satisfy 
all the required conditions of the continuity of the longitudinal 
displacements u(z,s) over the section z = constant. 

Each of the functions cp. (b) for a choice of the required 
magnitudes TJ^(z) thus has, by the method here indicated, a very simple 

geometrically evident expression; namely, it is different from zero only- 
on rectilinear parts of the contour meeting at the Joint i; within the 
limits of each of these segments it is represented as a function of s by 
a linear function of s and assumes a value equal to unity at a point 
coinciding with the ith Joint, becoming zero at the other end point of 
the given segment. On all the remaining parts of the contour line the 
function <p^(s) will be identically equal to zero (fig. 2). 

It is clear that the chosen method of constructing the functions qpp(s) 
for the assumed model is not the only one. For the required functions Uq(z) 
there may be taken any m independent magnitudes. To each aggregate 
of m independent magnitudes U^z ) will then correspond an aggregate 
of m linearly independent functions cp^(s) each of which will be 

continuous over the entire multiply connected contour and on an individ ual 
segment of the contour is represented by a linear graph. Thus, for 
example, if we wish in the first of the series (l.l) to separate initially 
the longitudinal displacements relative to the elementary computation of 
the shell as a beam of a composite multiply connected cross-section on 
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the "basis of the hypothesis of Bernoulli assumed for the entire section 
z = constant, we must, for three of the m possible displacements Uq(z) , 
take the values Uq(z ), Up ( z ) , and z ) determining the displacements 

of the hinged model on the prismatic surface as a rigid plane system. 

The functions cp^(s), cp£>(s), and cp^s) corresponding to these magnitudes 

are each linear functions of the cartesion coordinates x = x(s) and 
y = y(s) of a point of the contour of the cross-section of the shell. 

The remaining terms of the series (l.l) for the function u(z,s) determine 
a state of longitudinal displacements for which the cross-sections of 
the shell, after deformation, do not remain plane. Retaining the terminology 
adopted in our work on the general theory of thin-walled rods (reference l) 
we shall denote the deviation from the law of plane cross-sections as 
"deplanations" of cross-section. 

Hie deplanation of a multiply connected shell is thus determined by 
the m — 3 independent quantities U^(z), U^(z), ..., U^z). These 

magnitudes will be called the generalized coordinates of the deplanation. 

In constructing the functions ^(s) (k = 1,2, 3 , . .., n) entering 

equation (1.2) for the transverse tangential displacement v(z,s) of a 
point of the contour (z,s), we shall start from the deformation of an 
elementary transverse lamina (fig. l(a)) in its plane z = constant. 
Considering this lamina as a rod system and assuming the elements of the 
frame to be inextensible we arrive at the result that the contour displace- 
ment v(z,s) may be expressed in terms of the displacements Yjc(z ) 

(k = 1,2, 3 , n) of the plane kinematic model considered above in the 

plane of the cross— section of the shell. The coefficients of Y^z) of 
the sum (1.2 ) can be taken as independent magnitudes determining the form 
of the displacements of the hinged rod system in its plane. The number n 
of required functions Yjj-(z) is equal to the number of degrees of 
freedom of this system in the plane of the cross-section and is determined 
by the formula 


n = 2m — c ( 1 . 3 ) 

where m is the number of joints and c is the number of rods of the 
transverse multiply connected elementary frame. 

By choosing in some manner n Independent magnitudes Y^ ( z ) for 
the displacements of an elementrary rod system in the plane of the 
section z = constant and giving in succession to each of these magnitudes 
unit values and assuming the others equal to zero we can, by considering 
the elementary displacements thus obtained of this system, determine al 3 
the required functions "'t'k(s). Each of these functions will represent a 
contour displacement of a point s in correspondence with the elementary 
state Y k * = 1 and V^* = 0 for h ^ k (fig. 3). The function ty k (s) 

within the limits of each straight segment of the contour of the shell 
maintains a constant value (does not depend on s) and represents an 
axial displacement of the corresponding hinged rod model. We shall 
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thus set up n linearly independent diagrams of the functions tjj-(s) far 
any choice of the magnitudes Y z (z). 

These magnitudes will be chosen in such a manner that three of 
them, Yq ( z ) , Yp(z ) , and Y^(z ), will refer to the displacements of the 

model as a plane system of rods as a whole without change in shape of this 
system. The remaining magnitudes Y^z), Y^(z), ..., Y n (z) will refer to 
such displacements of the system for which the mutual position of the 
individual links of the system changes. The phenomena associated with 
the change in 3hape of the hinged rod system will be denoted as the 
deformation of the contour of the system. The contour deformation is 
thus determined by n — 3 independent magnitudes where n is the number 
of degrees of freedom of an elementary lamina of the shell considered as 
a plane hinge-connected Bystem. The magnitudes Y^(z ), Y^(z), ... Y n (z) 

will be denoted as the generalized coordinates of the deformation of the 
shell contour. 

The functions ljr-^s) chosen in some way or other corresponding to 

the n degrees of freedom of the system of rods in its plane satisfy the 
condition of linear independence and the condition of continuity of the 
displacements determined by these functions of the elementary transverse 
frame at all points of its contour including also the nodal points since 
the hinged model in each of the n possible elementary states Yv-* = 1 " 

Y] t (z) = 0 for h ^ k remains everywhere continuous. 


2. Fundamental Differential Equations — Method of Displacements 

After choosing the functions qpq(s) and ^(s) in the sums (l.l) 
and (l.2) the problem reduces to the determination of the functions TJ^z) 
and Y fc (z) (i = 1,2,3, . . . , m, k = 1,2,3, ..., n). 


Let cr = cr(z,s) and t = t(z^s) denote respectively the normal and 
tangential stresses arising in the section z = constant. We shall consider 
these stresses functions only of the coordinates of the point (z,s) on the 
surface assuming that over the thickness of the shell the stresses o 
anl t are uniformly distributed. On the basis of Hooke* s law 


a 



T 




Hence making use of equations (l.l) and (l.2) we have 


( 2 . 1 ) 


ct(z,s) = E ytzJcp^B) 


(i = 1,2,3, m) 


( 2 . 2 ) 


NACA TM 1234 


5 


t(zjs) = G 


Ui(z )qp ± *(s ) + y~ 



1,2,. . . , Htjk — 1,2,..., ll) 

( 2 . 3 ) 


where E and G are the moduli of elasticity of the material of the shell 
in tension, and shear. 

The elementary cross— sectional lamina of the shell will "be under the 
action of external forces consisting of normal and shearing forces acting 
in the sections z = constant and z + dz = constant and the given 
surface forces. Let p* = p*(z,s), and q* = q*(z,s) denote the forces 
external with respect to the given lamina acting respectively along the 
generator of the shell (positive in the direction of increasing coordinate z) 
and tangential to the contour line of the cross-section (positive in the 
direction of increasing coordinate s). Referring these forces to unit 
area of the middle surface we obtain the equations 

= |^8 + 1 ( 2 . 4 , 2 . 5 ) 

8 = 8(s) is assumed a given (in the 
general case not continuous) function of a single coordinate s and the 
magnitudes p = p(z,s), and q = q(z,s) represent the given external 
surface forces. 

The integral conditions of equilibrium of the elementary frame for 
chosen forms of the displacement determined by the m + n degrees of 
freedom, on the basis of the principle of Lagrange can be represented in 
the form of m + n equations: 


/If V* ~/ T<P J di ’ + / ptp 3 d8 = 0 (J = 1,2,3,..., rn) (2.6) 


-X3 /^ d8 + /* h ds = 0 1 , 2 , 3 ,..., n) ( 2 . 7 ) 


p* = ~8 + p 

where the thickness of the shell 


where dE is the differential of the area of the shell cross-section 


dF = 8ds 


( 2 . 8 ) 
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The integrals above and those given in the following are definite and 
are taken over the entire contour of the section z = constant. 

Equations (2.6) represent m conditions of equilibruim of the 
lamina dz = 1 in the direction perpendicular to the plane z = constant; 
equations (2.7) represent n conditions of equilibrium of the same 
lamina in the plane z = constant. Each of the equations (2.6) expresses 
the equating to zero of the sum of the work of all external and internal 
forces of an elementary lamina with change in the deformed state of the 
lamina relative to its plane. 

For the virtual displacements in the equations of subscript J of 
this group there are taken the longitudinal displacements uj 

the points of the elementary frame determined only by the term of 
•subscript j of the sum (l.l) for Uj* = 1. The first and third terms 

of this equation refer to the work of the external forces p* = ^ 5 + p 

acting on the element of width dz = 1 and directed perpendicular to 
its plane. 

The middle term expresses the work of the internal shearing forces. 
For an element ds this work is determined as the product (with reversed 
sign) of the shearing force tS ds by the shear deformation which in the 
case considered of the variation of the deformed state is equal to the 
derivative cpj*(s) of the function cp^(s). 

Each of equations (2.7) was obtained by equating to zero the sum of 
the work of all the external and internal forces of the elementary frame 
on the corresponding displacements for a change in the deformed state of 
the element in its plane. For the virtual displacements in the equation 
of subscript h there were taken the transverse contour displacements 
Vh = 1 r h( 3 ) of the element determined only by the term of subscript h 
of the sum (1.2) for the generalized coordinate Y^* = 1. The first 
and third terms of equations (2.7) refer to the work of the external 

St 

contour forces q* = ^ 5 + q of the element acting in its plane. The 

second term expresses the work of the internal forces on the deformations 
of the element corresponding to the h— th elementary state of the 
displacements of the hinged system in its plane. For the element ds 
this work in the case of bending is determined (with reversed sign) as 
the product of the bending moment 



M(z,s) = \( z ) M k;( s ) 

Mu(s) 

by the angle of rotation — ds 

lUJ 

element. The letters = M^( s ) 


(k = 1 , 2 , 3 } n) (2.9) 

of two adjacent sections bounding this 
and = Mj-^s) denote the bending 
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moments of the transverse frame element corresponding to the elementary 
states of deformation of this frame Yk* = 1 and Y^* = 1. These moments 

are found by the usual methods of structural mechanics by proceeding 
-backward from the deformed state of the rod system to the internal forces. 

The magnitude J = j(s) represents the moment of inertia of an 
arbitrary cross-section of the elementary frame of the shell having a 
•width dz = 1. If the shell consists of only rectangular plates and has 
no transverse connections then evidently 

J = ^ (2.10) 


where S is the thickness of the corresponding plate. 

In the case of a shell reenforced with additional transverse frames 
the magnitude J must be computed with account taken of the mean moment 
of inertia of these frames, that is the moment of inertia associated 
with a unit length of the shell. 

Substituting in equations (2.6) and ( 2 . 7 ) for a and t their 
expressions ( 2 . 2 ) and ( 2 . 3 ) we obtain a system of m + n linear differential 
equations with respect to the required generalized displacements, m 
longitudinal Uq(z) (i = 1 , 2 ,..., m) and n transverse Yjj.(z) 

(k = 1,2,..., n). This system can be represented in the form 


7 21 a Ji u i" “ 21 b ji u i “ 21 c jkV + g P J = ° 
i i k 

/i,J = 1 , 2 ,..., m\ 
\h,k = 1 , 2 ,..., n/ 

21 c mV + 2T r hkV - 7 21 BuJk + k% = 0 

i k k 


( 2 . 11 ) 


where 


7 is a constant magnitude defined by the equation 


7 = 


E 

G 


( 2 . 12 ) 
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The coefficients of equations (2.1l) are computed by the formulas: 


a ji = J ‘PjUtyi c jk = J qp j 1 ( s H-jj- ( s ) aP, r hk = / ^h^ s ^k^ 8 ^ 35 


(2.13) 


bji = J <p J *(s)cp 1 , (s) dP, c M = J ^(sJqPi^s) dP, s 


'Mh(s)%(s) 


hi: 


EJ 


da 


where the integrals are taken over all the elements of the cross-section 
of the shell. These coefficients possess the properties of symmetry 


a ji “ a iJ* 
*Ji = 


^hk = r kh 

„ c hi 

B hk 8 kh 


c^jj. for h = k 


(2.14) 


expressed by the theorem of Betti on the reciprocity of the work of an 
elastic system. 

Formulas (2.13) are of a general character and permit computing the 
coefficients of the equations (2.1l) for a shell of arbitrary contour of 
cross-section for any method of approximation of the required displace- 
ments u(z,s), t(z,b) for -varying s. 

In choosing the functions 9^(3) (i = 1,2,..., m), ^(s) (4 = 1,2,..., n 

by the method described above the quadratures on the right sides of 
formulas (2.13) for each straight segment of the contour receive a simple 
expression because each of the functions cpj_ over this segment depends 
on the coordinate s linearly and the derivative <p like each of the 
functions has a constant value over a straight ' part of the contour. 

The quadratures of the first five formulas (2.13) for dF = S ds 
have the same form as the last of these formulas except that instead of 
the magnitudes reciprocal to the moment of inertia the thickness of the 
shell enters under the integral sign. All the coefficients of the equations 
can be ,computed by the known devices of the theory of framed structures 
with the aid of diagrams of the functions 9^(3), 9 i *(s), ^(s), M^s) 
constructed for the entire multiply connected contour. 

The magnitudes p j(s) and qj 1 (z ) referring to the free terms of 

equations (2.11 ) are known functions of z and for given surface forces 
of the shell p(z,s), q(z,s) are computed by the formulas 
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pcpj da, q^ = I q%, ds (2.15) 

where the contour integrals are taken over the parts of the contour for 
•which the expressions under the integral sign are different from zero. 



The magnitudes pj ( z ) and 

meaning arising from tjie method of their definition can "be termed the 
unit (referring to unit length of the shell) generalized external forces. 
The longitudinal one, p^(z), is computed as the work of the external 
longitudinal surface forces, p(z,s), over the longitudinal displacements, 
<Pj(s), of the elementary (unit) state of deformation of the lamina 

for Uj* = 1, while the transverse one, q-^(z), is computed as the work of 

the external surface contour forces q(z,s) obtained for the contour 
displacements of the points of the elementary lamina 4 r } 1 ( s ), determined 
by the unit displacement U^* =1. 


q^( z ) corresponding to their physical 


Formulas ( 2 . 13 ) and ( 2 . 15 ) are easily extended also to shells 
stressed over individual longitudinal elements (stringers) and under the 
action of concentrated loads in the section z = constant. In this case 
the quadratures in (2.13) and ( 2 . 15 ) must be taken in the sense of the 
Stieltjes integrals. To the integrals for continuous distribution over 
the contour of the differentials dF, p ds, q ds must be added magnitudes 
which represent the sum of the products of the finite factors concentrated 
at definite points (the areas of the stringers, the concentrated forces) 
and the values of the corresponding functions under the integrals at 
these points. 


The systems of the fundamental differential equations (2.1l) of the 
prismatic shell may be represented in the form of table 1 ( 2 . 16 ) (found at 
the end of the text). This table represents the differential matrix and 
the free terms of equations (2.1l). The elements of the matrix represent 
the linear differential operators with constant coefficients, where these 
operators in the first and fourth quadrants are of the second order and in 
the second and third quadrants are of the first order (D and Jp denote 
the derivatives with respect to z of the first and second order of the 
corresponding required function written in the top row over the column 
containing the given differential element of the matrix). To obtain 
any differential equation of the system it is necessary to multiply each 
of the elements of the corresponding row of the matrix symbolically (by 
which is meant the differential operation expressed by the given element) 
‘by the function in the top row over the column containing the given 
element, combine all these products and to the sum add the free term of 
the given row and equate the result to zero. The system of differential 
equations represented in the form of table 1 consists of m + n 
ordinary linear differential equations with constant coefficients each 
of the second order with respect to the m + n required functions 
U* (z ) and V v (z). 


I I HI ■■■II II ■ ■ ■■” 
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This system of equations by its physical meaning divides- itself 
into two groups of equations. The first group written out in the first 
half of the table consists of m equations and expresses the equilibrium 
of the elementary elastic frame lamina out of its plane as a system of 
rods possessing on the prismatic surface m degrees of freedom. The 
second group of equations, consisting of n equations, refers to the 
equilibrium of this frame in its plane possessing in the latter n degrees 
of freedom.. The matrix of the differential equations possesses a 
symmetrical structure and consists of four quadrants of which the first 
and fourth are the principal quadrants, the second and third the secondary, 
and represent energostatically the mutual effect of the two forms of 
deformations of the elastic shell. The elements of the secondary 
quadrants of the matrix which are symmetrical with respect to the 
principal (diagonal) terms differ only in sign. The equality of these 
elements in absolute value is a consequence of the theorem on the 
reciprocal work of elastic systems. 

The equations here described represent a generalization of the 
earlier obtained eight— term equations for cylindrical shells and 
composite systems of open profile (references 1 and 2). The equations 
represented in the form of table 1 are completely analogous to the 
equations of the theory of statically indeterminate systems and differ 
from the latter in that these equations are differential since the 
deformed state of the transverse plane frame as an element of the shell 
depends on the position of the element along the length of the shell 
(on the coordinate z). 

The differential equations (2.1l) were derived for an arbitrary 
choice of the functions cp^s) and ^(s) dete rminin g according 

to equation (2.13) the coefficients of these equations. Since the functions 
qPj- and are linearly independent and each of them may be given with 

an accuracy up to an arbitrary factor, for the required functions TL^z) 
and V k (z) there can always be chosen such independent generalized 

longitudinal and transverse displacements of the elementary hinged 
system of the shell dz = 1 for which the polygonal functions ^(s) 
and ^ k (s) corresponding to these displacements over the entire cross- 

section z = constant possess the property of orthogonality. Assuming 
each of the pair of functions cp^s) and i|r k (s) orthogo nal , we hVihIT 
have 


a Ji ~ 


cpjqpi dF = 0, 


if j/i 


rhk = 


♦h+k 31 = °> 


if 


h f k 


( 2 . 17 ) 



HACA TM 1234 


11 


For these conditions equations (2.1l) assume the form 


7a JJ U j" “ h 21 U ± ~ ^ c JkV + h Pj " 0 = i* 2 *---* m ) 

i k 

( 2 . 18 ) 

XT c hi u i* + r hh T h" - ?II V T k + g <3h = 0 = 1,2,..., n) 

i • k: 


The differential matrix of these equations now has a simpler 
structure since, with the condition of orthogonality ( 2 . 17 ) satisfied, 
the secondary differential terms of the principal (first and fourth) 
quadrants, table 1 , that is the terms with the coefficients a^j and rj^, 

for i ^ j, h ^ k vanish. 


Written out in full equations (2.18) are represented in table 2 (2.19) 
(found at the end of the text). The magnitudes U^(z ) and V^z) corre- 
sponding to the conditions of orthogonality may be called the principal 
generalized coordinates of the longitudinal and transverse displacements of 
the shell. 

The choice of the orthogonal functions 9 ^( 0 ) and ^(s) may be 

made by the graphico— analytical methods of structural mechanics by 
drawing the diagrams of these functions and orthogonalizing the group 
states of the elementrary displacements of the transverse laminar element. 


3 . Boundary Effect — The Internal Generalized Forces — 

Longitudinal and Transverse Bimoments 

The equations (2.19) or, in the case of an arbitrary choice of 
the functions 9 ^ and t)^, equation ( 2 . 16 ) possess the same structure 

as the equations of vibration of an elastic system possessing a finite 
number of degrees of freedom. Of the present day methods of integration 
of a symmetrical system of linear differential equations with constant 
coefficients the most effective method is that of A. N. Krylov, which 
permits reducing this system rapidly to its equivalent single differential 
equation. In our case this equation will be of the order 2(m + n). It 
follows that the required functions Ej^(z) and Vj c (z ) satisfying the 

system of equations ( 2 . 18 ) will be determined with an accuracy up to 
2(m + n) arbitrary constants. The number of these constants is equal to 
twice the number of degrees of freedom of the elementary lamina of the 
shell dz = 1 in space. This is in full agreement with the number of 


12 


NACA TM 1234 


independent kinematic conditions which can he given for the end sections 
of the shell z = 0, z = l (where l is the length of the shell in the 
direction of the generator). The position of all the points of any of 
these sections in space for a shell treated in the section z = constant 
as a discrete system is determined by m + n independent magnitudes of 
which m magnitudes U^fz) determine the position of 1 points of 
Juncture of the section on a prismatic surface (along the generator of 
the shell) -and n magnitudes V k (z) determine the position of these 

points in the plane of the cross— section. For a single end section of 
the shell m + n arbitrary magnitudes may thus be assigned. For the 
two sections z = 0 and z = l which bound the given shell over its 
length the number of independent conditions is equal to 2(m + n) which 
corresponds to the number of arbitrary constants in the integrated 
equations (2.11 ). 

With these constants given we can obtain a solution for a given 
shell for the most varied boundary conditions relative to the longitudinal 
and transversed displacements, this solution being entirely determinate 
and unique and satisfying all the required kinematic conditions of the 
shell . 

Let us now consider the problem of the equilibrium of the shell for 
which the boundary conditions in the sections z = 0 and z = l are 
given in terms of stresses or, in the case of the mixed boundary problem, 
partly in terms of stresses and partly in terms of displacements. By 
definition of the functions IL^(z) and Y-^( z ) from equations (2.1l) 

the stresses a = cr(z,s) and t *= t(z,s) at any point of the cross- 
sections z = constant will be found according to equations (2.2) and 
(2.3) also with an accuracy up to 2(m + n) arbitrary constants. The 
stresses a and t at the section z = constant for chosen polygonal 
functions cpq(e) and ^(s) can be expressed in terms of m + n 

independent generalized static magnitudes. These magnitudes we shall 
introduce by generalizing the fundamental concepts of the elementary 
theorem of the bending of beams aB was done in our work on the general 
theory of the thin-^walled elastic rods based on the law of sectorial 
areas (reference 3). 

Starting from the idea of the virtural work of the normal and 
shearing forces 06 and t 5 of the cross-section z = constant on 
each of the m + n possible displacements of the points of thiB section 
in space we introduce the following magnitudes 



fcxcpj dF 

(j 1^2,..., in) 


«h ( s > -J 

1 Ti|f h dF 

(h = 1,2,..., n) 

(3.1) 
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where the contour Integrals are taken over the entire cross-sectional 
element of the shell. The magnitudes Pj, represent generalized 

longitudinal and transverse forces of the section z = constant. 
Considering these magnitudes as the internal forces of the shell we 
express them in terms of the fundamental functions and Y^. On 

the "basis of equations (2.2), ( 2 . 3 ), and (3-l) we have ~ 




= E 


^h = 


= Gr 


a jiV 


C hi U i + 


(£;5 1 0 (3 - 2) 


3 vhkTk* 


If the functions cp^ and ^ are orthogonal 


P J - 


Qh = G (^~ + r hhV^ (j — m) (3«3) 


Equations (2.2) can now "be written as 


or in expanded form 

( \ P l( z ) t \ ?2( z ) r \ 

a(z,s) = -J* qh(s) + — cpo(s) + 

» a -, -[ aoo e ~ 


a22 


,, m) 


p m( z ) 

• + <p m ( s ) 


a Tnm 


(3.^) 


(3.5) 


Equation (3.4) represents the generalization of the three— term 
formula proposed "by us for the normal, stresses of open cylindrical wnfl 
prismatic shells (for thin-walled rods with open profile sections) 
undergoing simultaneous tension, "bending in two planes and torsion. The 
theory of such shells is based on the hypothesis of the nondeformability 
of the contour of the cross-section. 

« 

For the open shell with rigid cross-sectional contour setting 


cpi = 1, c p 2 = x(s), cp 3 = y(s), = o>(s) (3.6) 

where x(b) and y(s) are the cartesian coordinates of an arbitrary 
point of the contour in the principal central axes of the section 
and <n(s) is twice the area of the sector determined by the arc MqM 
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and the straight lines joining the ends of this arc with the center of 
flexure (fig. 4) we obtain for the generalized forces Pi,..., P^, the 

first four terms of the series (3.4), the following values (reference 3, 
p. 48). 


?1 = 


al dF 





P 


3 


ay dF 


Mr, 1*4 



cxu dF = B 


(3.7) 


The first three equations of (3*7) determine the known statical 
magnitudes (normal force and moments) of the cross-section of a beam; 
the fourth equation determines a new static magnitude having the 
dimensions kilograms per centimeter^ and representing the work of all the 
elementary longitudinal forces a dF of the shell on the deplanation of the 
cross-section determined by the law of sectorial areas. This magnitude is 
termed by us the bimoment. The condition of orthogonality of the four 
fundamental functions ( 3 * 8 ) is determined by us in the form 


lx dF = / ly dF = / xy dF = 0, lto dF = / xo) .dF = / yiu dF = 0 ( 3 . 8 ) 


Of these conditions the first three coincide with those known from 
the theory of the bending of beams and which determine the principal 
central axes of the cross-section of the shell. The second group of 
equations ( 3 * 8 ) refers to the sectorial geometric characteristics of 
the cross— section of the shell and determines the sectorial origin M Q 
(the starting point for computing the sectorial area) and the center of 
flexure of the shell. For the conditions ( 3 . 8 ) the bimoment B = B(z) 
represents a generalization statically equivalent to the vanishing of 
the longitudinal force. 

Corresponding to equations (3*6) and (3-8) we obtain for the 
geometric characteristics ajj (j = 1 , 2 , 3,4) of the first four terms 

of the series ( 3 - 5 ) the values 



(3.9) 


The first three of these characteristics agree with the well-known 
fundamental characteristics of the cross-section (area and moment of 
inertia) in the theory of the bending of a beam. The fourth charac— g 
teristic represents a new magn itude having the dimensions of centimeters 
and denoted by us, in analogogy with J x , J , the sectorial moment 
of inertia. 
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The equation (3*5) for the thin-walled open shell, keeping only the 
first four terms, assumes on the basis of equations ( 3 • T) a^-i (3*9) the 
form 




F 


x + 


7 


B 


0) 


j cd 


(3.10) 


The first three terms of this equation determine the stresses of 
the thin-walled rod in the case of tension (compression) and bending. The 
fourth term determines the stress arising in the case of torsion of the 
rod and distributed over the section by the law of sectorial areaB <jd = cd(s). 

z), and = 

by the well-known methods of the resistance of materials. The static 
magnitude B = B(z) refers to the normal stress = (B/j^)^ due to 

the torsion determined from the differential equation 


The static magnitudes N = N(z), Mj. = M^ 


l^.(z) are found 


B" - 




EJ, 


B + 


m 


CD 


EJ ( 


= 0 


CD 


(3.U) 


in which GJ^ is the stiffness of the rod in pure torsion determined by 

the theory of Saint Tenant and m = m(z) is the external torsional 
moment relative to the center of flexure. 

From the analogy here given it follows that equation (3*5) for the 
shell with closed multiply connected deformable contour for a choice of 
thb first three functions <Pj(b) by the equations cp^ = 1, = x(s) 

and cp^ = y(s) and the conditions of orthogonality (2.1 f) can be 
represented in the form 


M 




c - 1 + 5t 1 + r y 


7 


a kk * ^5 5 


P 

m 

+ cp 

a. L ... m 
vmm 


(3-12) 


where the first three terms refer to the stresses distributed over the 
section according to the law of plane deformation and the corresponding 
elementary theory of the bending of beams. The remaining generalized 
longitudinal forces P] 4 .(z), P 5 (z), ..., P m (z) for the contour of the shell 

possessing in the longitudinal direction a number of degrees of freedom 
greater than three will represent the internal "longitudinal" forces 
having essentially the same character as the forces N, M^, and My- of a 

beam, the only difference being that these forces, for a choice of an 
orthogonal system of functions <p^(s), each represent over the cross- 
section a balanced system of longitudinal forces and arise as a result of 
the deplanation of the section. In contrast to the longitudinal force N 
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and the tending moments and My. the generalized longitudinal 

forces P^ , . . . , P m associated with the deplanation of the section 

will he denoted as the longitudinal bimoments. These bimoments 
correspond to the generalized coordinates (f)^, cp m of the 

deplanation of the cross-section. 

The geometric characteristics 



by analogy with the well-known magnitudes 


( 3 . 13 ) 


a22 = J x 2 dP = Jj, a^3 = J y 2 dP = J x ( 3 . 1 ^) 

are denoted as the longitudinal bimoments of inertia. 

In a similar manner the physical sense of the generalized trans-verse 
forces Q^(z) (h = l, 2 ,...n) determined by the corresponding equations ( 3 .l) 
and (3.3) can be explained. 

If of the n generalized transverse displacements we give the 
first three functions Yj(z)j Vg(z), ^ V^z) the sense of three 

independent displacements of the elementary lamina as a rigid system in 
the plane of the cross-section of the shell and set 

^(b) = x* ( s ) j t 2 (s) = ^3(3) = x(s)y*(s) - y(s)x*(s) (3-15) 

the magnitudes Qq_(z), Qg(z ) , and Q^( z ) corresponding to the coordinates 

of the displacements will represent the first two transverse forces and 
the torsional moment of the entire section z = constant , respectively. 

The remaining static mag n itudes Qij.(z), Q^(z), and ^(z) will represent 
the generalized transverse forces each computed as the work of the 
elementary shearing forces t dF by the corresponding contour (tangential) 
displacements j that is by the displacements of the shell determined by 
the generalized coordinates ^(s), ^(s),...^ \jr n (s) of the deformations 

of the contour of the cross-section. 

In contrast to the forces and Qg = Qy and the torsional 

moment Q3 = H, the generalized forces corresponding 

to the components Vlp v 5 , • • . , Y n of the deformations of the contour of 
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the section we shall denote as the transverse hi moments of the shell. 
Equations (3*3) determine the relation between the generalized forces 
and the generalized displacements of the shell. 

The generalized forces Pj(z) (J = l* 2 ,..., m), Q^(z) (h = 1,2,... , n) 

characterizing the state of the normal and shearing forces in the cross- 
section of the shell are likewise determined with an accuracy up to 2(m +, n) 
constants of integration of the system of differential equations of the 
shell. We now assume that any edge of the shell z = z Q is acted upon 
by a given system of normal and shearing forces. Let p°(z 03 s), q°(z 0 ,s) 
represent the corresponding normal and shearing forces referred to unit 
length of the contour at the point s. These forces as functions of 3 
may be given entirely arbitrarily. In our problem the longitudinal and 
transfer se forces at any section z = constant are determined by a 
finite number of independent static magnitudes of the corresponding 
finite number of degrees of freedom of the displacements of the contour 
line around the section and in its plane. These magnitudes are gener- 
alized longitudinal and transverse forces which for the given normal and 
shearing forces p°(z 0 ,s), q°(z 0 ,s) are determined by the equations 


Pj° = /p° 9 j < 1 b, Qh° = f q% ds (3 = 1 , 2 ,..., 

The above equations together with equations (3.2) 
section z = z Q , to the equations 

Ea Ji U i f = J P 0c Pj ds > G ( c hi u i + 


m, h = 1 , 2 , . . ., 
lead, for the 



ds 


n) (3.16) 


(3.17) 


These equations determine the relation between the required generalized 
displacements and the given generalized forces at the bounding 
section z = z Q . 

Having the general integral of the differential equations of the 
shell and making use of the general formulas ( 3 . 17 ) we can determine 
the state of the deformations and stresses of the shell for the most 
varied boundary conditions on the sections z = 0 and z = l given in 
terms of stresses, displacements, or partly stresses and partly displacements. 


b . Shells Having a Single Degree of Freedom for Deplanation 

of the Section and Deformation of the Contour 

Let us consider the class of shells for which the deplanation of 
the section and the deformation of the contour are each determined by a 
single parameter which is a function of the coordinate z. To such 
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shells may he reduced a number of practically important problems encountered 
in the construction of thin-^walled structures in various applied fields of 
technology. 

As a very simple example let us consider the shelly the cross-section 
of which is represented in figure 5* We assume that the vertical plates 
of the shell have cylindrical hinges at the lover longitudinal edges and 
are fixed over the entire supporting length against longitudinal displace- 
ments. For such a shell the elementary transverse lamina as a plane 
hinged system with the supporting points completely fixed, possesses 
three degrees of freedom, two out of its plane and one in its plane. 

We assume further that the section of the shell is symmetrical relative 
to the middle vertical line and that the shell is acted upon by a 
transverse load of given intensity q = q(z ) applied in the plane of 
the horizontal plate. In this case the longitudinal displacements of 
the shell u(z,s) in the section s = constant for any two points of 
the section symmetrically located with respect to the axis of symmetry 
will be equal in absolute value and opposite in sign. 

These displacements are determined by a single parameter Up(z) 
for which we take the longitudinal displacement (taking it positive) of 
the upper right angle of the section. The graph of the displacements 
represented for Up* =1 as the function q>^(s) is shown in figure 6. 

The law of distribution of the longitudinal displacements and therefore 
also the normal stresses over the section z = constant is postulated 
by the displacements of the points determined by this function. This 
law differs from that of plane sections. 

The derivative cpp* (s ) within the limits of each of the segments 

of the section remains constant and in absolute value equal to 1 /dp 
for the vertical portion and 2/cLp for the horizontal portion where 
dp and dp are the lengths of the vertical and horizontal portions. 

The transverse contour displacements of the shell for the case 
where the longitudinal edges are fixed, are determined likewise by a 
single parameter Vp(z). For this parameter we take the bending of the 
horizontal plate in the section z = constant (fig. 7)» The contour 
axial displacement of the points of the section of the shell that are 
determined by the function i^p(s) corresponding to unit value of the 

required magnitude Tp* =1 is shown in figure 8. These displacements 

on the vertical portions of the section are equal to zero and on the 
horizontal equal to unity. The arrow in figure 8 gives the direction 
of the positive displacement. Figure 9 shows the graph of the bending 
moments of the elementary transverse frame corresponding to the 
displacement Tp =1 of the horizontal elements. 

For the moment Mp at the joint of the frame it is not difficult 
by the methods of structural mechanics to obtain the formula 
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2dn 


6 e 

J 1 + d l 


^l' 


(4.1) 


where dq and dg are the widths of the vertical and horizontal plates, 
and J 2 , the vertical and horizontal moments of inertia of the 
elementary frame. 


In the problem under consideration there are thus to he determined 
two fundamental functions, the longitudinal displacement Uq(z) 
determining the deplanation of the section and the transverse displace- 
ment Vq(z ) determining the deformation of the contour of the section. 


Setting in equations (2.1l) i, J = 1; h,k = 1 we obtain 


7 a n u i" ~ ^n 0 ! ~ c n v i* = c n u i* + r u T i" ~ 7 a ii v i + | ii 


0 (4.2) 


The coefficients of equation (4.2) are computed by formulas (2.13). 
Applying these formulas and making UBe of the graphs of the corresponding 
functions given in figures 6 , 8, and 9 we obtain 


*11 “ 




dF = | (2 F-l + F 2 ), 




r ll 



dF = F 2 


(4.3) 


'Mq 2 ds 


=11 = 



where 5^ and S 2 are the thickness of the vertical and horizontal 
plates respectively, Fq and F 2 , the cross-section areas of these 
plates and and J 2 , the moments of inertia referred to unit 

length of the longitudinal sections of the shell on the vertical and 
horizontal parts, respectively: 


Fq ^ dqSq, 


Sq3 

Jl = IP 


F 2 “ d 2 5 2 
So 3 


Jo = 


12 


(4.4) 

(4.5) 
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If the shell has along its length a number of transverse frames 
the mome nts of inertia Jq and must "be computed with account 

taken of the resistance of these frames to tending. 

Equations (4.2) for the coefficients (4.3) will be completely 
determinate for a given shell. By giving different values to the 
magn itudes dq, d^, Sq, Sq’ Jq* J 2 a number of special problems 

can be considered referring to the computation of a given shell for 
an antisymmetric load on the cross-section. Leaving in equations (4.2) 
the coefficients a ^ ) b-^,..., Sqq as yet arbitrary let us consider 
more in detail the solution of these equations and incidentally explain 
the physical side of the elastic computation model. 

The system of differential equations (4.2) by eliminating the 
longitudinal displacement reduces to a single equation relative to the 
transverse displacement of bending of the horizontal plate Yq(z). This 
equation may be written in the form 


where 



- 2A 2 Y 1 " 



+ - ^ a ni’’ ) = 0 ■ 

c n G 


(4.6) 


A 2 - 


4 a n r n 


^ a n s n + B ii r ii - °u 2 ^ 7 = 1 (4,7) 

a n r n G 


Differential equation (4.6) expresses the equilibrium of the elementary 
transverse frame in the dense elastic medium which is constituted by the 
given shell relative to this frame. In its form this equation agrees 
with the equation of the bending of a beam on an elastic support resisting 
not only bending but also shear along the line of contact of the beam 
with the support. The elastic medium in equation (4.6) Is represented by 
two independent magnitudes A 2 and B , which may be called the generalized 
elastic characteristics of the beam, each having a single degree of 
freedom for the deplanation of the section and the deformation of the 
contour. These characteristics, as is seen from equations (4.7) and. (4.3) 
are determined by the functions cp^( s ) and Y-^(s) referring respectively 


to the deplanation and deformation of the contour of the shell and by 
the form and geometric dimensions of the cross-section of the shell, by 

Tj 1 

the ratio 7 = % of the elasticity modulus and the stiffness in bending 

in the transverse direction represented by the moments of inertia Jq 
and Jp of the elementary frame of width dz = 1. Thus, for example, 
for the shell represented in figure 5 the elastic characteristics A^ 
and for given dimensions of this shell on the basis of formulas (4.3) 

assume entirely definite values. 
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The characteristic B^ aa is seen from the second of formulas (4.7) 
is proportional to the coefficient B y] referring to the tending of 

the elementary transverse lamina as a frame. If the plates of the shell 
are^hinge connected at the joints, the transverse tending moments of the 
elementary lamina on changing its shape is equal to zero. The coeffi- 
cient s -| -| and therefore also the characteristic b 4 in this case 
Becomes zero: 


8 11 = 


B^ 


0 


(4.8) 


and equation (4.6) for the zero moment shell goes over into the 
simpler equation 


Vi IV _ SA 2 ^" + ^gdnq - ra-^q") = 0 

in which now 

o ^l^n “ c u 2 

— - - 

4a ll r ll 


(4.9) 


(4.10) 


Proceeding to the solution of the differential equations of the 
class of shells here considered we note that the two kinematic magni- 
tudes Bq(z) and Yq(z) ty formulas (3*2) are reduced to two kinematic 

magnitudes corresponding to the static magnitudes P^(z) and Qq(z) 

introduced ty us, the first of which represents the required longitudinal 
force generalized for the deplanation of the section (longitudinal 
timoment) and the second the generalized transverse force in deformation 
of the contour (the transverse timoment). We thus have four required 
magnitudes Uq(z), Y 1 (z), Pq ( z ) , and Qq(z ) representing generalizations 

of four magnitudes that are well known in the theory of the resistance 
of materials on the tending of teams: the deflection y(z), the angle 

of rotation of the section y*(z), the tending moment M(z), and the 
transverse force Q(z). 

The artitary constants of the integration are determined ty the 
Boundary conditions which, depending on the type of problem may Be 
given in generalized displacements Uq and Yq, in generalized forces Pq 

and Qq, or partly in displacements and partly in forces, only two condi- 
tions Being required at each of the Bounding sections of the shell z = 0, 
z = l. Thus for example if the infinitely long shell with section 
represented in figure 10 is acted upon ty a concentrated transverse 
horizontal force P the conditions for dete rmining the constants of 
integration on the part 0 <; z “ assume the form 



22 


NACA TM 1234 


G(cnUi + ruTi*) = — ^ for z = 0 

(4.11) 

G( c hUi + r-, iY-, 1 ) = 0 for z = » 


The first and third conditions are kinematic and express the fact that 
in the plane of action of the force P the deplanation of the section 
is equal to zero (section z = 0 by symmetry remains plane) and at 
infinity the deformation of the contour becomes zero. The second and 
fourth conditions are static and refer to the generalized (in the sense 
of virtual work) transverse force Q (the transverse bimoment) which 
in the section z = 0 must be in equilibrium, with the external transverse 
load and at infinity becomes zero. 

In table 3 are given the values of the generalized displacements 
for various sections of the infinitely long shell, represented on figure 10 
and loaded by a transverse horizontal force P. These displacements 
and stresses are computed for shells consisting of three plates of 
equal width and thickness for a ratio of thickness to width of plate 
equal to 0.01. 

In the first row of table 3 are given the values of the relative 

coordinate 5 = t (d is the width of a plate), in the second row are 
d 

given the ■ magnitudes proportional to the longitudinal displacements 
of the right joint of the "frame" of the cross— section; the third row 
refers to the transverse displacements Vq (deflection of the horizontal 

plate ), dete rminin g the deformation of the contour of the section of the 
shell. From these displacements are determined the transverse bending 
moments of the shell arising from the stiffness of the joint of the 
plates (horizontal and vertical). The maximum bending moments, according 
to the graph shown in figure 9 , arise at the joints of the transverse frame 
and are determined by the formula 

M(?) 

d 

in which the plus sign refers to the left and the minus sign to the 
right Joint. In the fourth row of the table are given the magnitudes 
proportional to the longitudinal normal stresses 0 = ff(£) referring 
to the points of juncture of the cross-section. The graph of these 
stresses over the section z = constant agrees with the graph of the 
deplanation of the section, shown in figure 6. The last row refers to 
the tangential, stresses t = t(£) forming in the section z = constant 
the flow of shearing forces s = ts. 
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The values of the displacements and stresses given in table 3 refer 
to the sections of the shell with positive abscissas £ = ^ = constant. 

For sections with negative abscissas the magnitudes Yq and a retain 
thefr" values while Uq and t reverse their signs. It is seen from 
table 3 that the deflections Y^ of the horizontal plate and therefore 

the transverse bending moments proportional to the deflections (the 
graph of these moments as a function of £ is shown in fig. 9) decrease 
at a slower rate with increase in distance from the plane of action of 
the force than the stresses a. 


5. Prismatic Shell with Closed Rectangular Section 

As a second example we consider a shell having in the cross-section 
a closed rectangular contour with two axes of symmetry (fig. ll). 

According to the theory developed above , the elementary transverse 
lamina dz = 1 of such a shell possesses 4° of freedom both in the plane 
and about the plane of the cross-section. The deplanation of the section 
and the deformation of the contour are each represented by a single parameter. 

Referring the cross-section of the shell to the principal central 
axes Ox, Oy the required displacements u(z,s) of an arbitrary point 
(z,a) of the middle surface is represented in the form 


u - Uqcp^ + XJgcpg + U3CP3 + TfyCfy, v = Y^Yq + ^2^2 + ^"3^3 + ^1^4 (5*1 ) 

where 


q>L = 1, q?2 = x(s), 93 = y(s), % = x(s)y(s) 

(5.2) 

Yq = h(s), i|r 2 = x*(s), Y 3 = y*(s), tl,. = x*(s)y(s) + x(s)y*(s) 

where x(s), y(s) are the coordinates of an arbitrary point of the contour 
determined by the distance s along the contour, h(s) is the length of 
the perpendicular dropped from the center of the rectangle on the line of 
the contour ; the primes denote derivatives with respect to the parameter s. 

The magnitudes U^(z)(i = 1,2, 3*4) refer to the generalized 
longitudinal displacements of the Bhell and for chosen functions cpq(s) 

(i = 1,2, 3, 4) the magnitudes Uq(z) are the displacements of the 

section z = constant in the direction of the Oz axis, the 
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magnitudes U 2 ( z ) and U^Cz ) are the angles of rotation of the section 

z = constant relative to the axes 0 y, Ox, respectively, and the 
magnitude U4 is the generalized deplanation of the section. The 
magnitudes V^Cz) (k = 1,2, 3,^-) represent: Y-^(z) — the angle of 

rotation of the section relative to the axis Oz (angle of torsion), 
Y 2 (z), Vo(z) — the deflections in the direction of the axes Ox 
and Oy, and V^(z ) — the generalized deformation of the contour of the 
section. 

It is easily shown that for chosen functions cpq(s) and ^(s) 
the system of differential equations (2.16) breaks down into four 
independent systems of equations. Of these systems three refer to 
tension (compression) and bending of the shell as a beam following the 
hypothesis of plane sections. This hypothesis in the first of 
formulae (5*1) is represented by the sum of the first three terms. The 
fourth system of equations consists of three simultaneous differential 
equations in the functions U^(z), Y^(z), and Y^(z). Leaving out in 

our problem the case of tension and bending we consider the equilibrium - 
of the shell in the case where it is acted on by external forces giving 
rise to torsion, deplanation of the section and deformation of the 
contour (fig. 12 (b)). 

In figure 13 is given the diagram of the function 94(e) = x(s)y(s) 
characterizing the law of variation of the longitudinal displacements 
for a unit deplanation Y4* = 1. This function represents a generaliza- 
tion of the law of sectorial areas previously given by us for thin- 
walled open rods. Figures lk and 15 give the diagrams of the 
functions = cp^*(s), Y-|_ = h(s). The signs of the magnitudes i|/^ 

and Y4 on these diagrams are denoted by arrows. The letters denote 
the length of the sides of the rectangle (fig. 11 ). Figure 1 6 gives the 
diagram of the transverse bending moments corresponding to the generalized 
deformation of the contour, taken equal to one. This diagram as well as 
the diagram of unit deplanation shown in figure 13 is antisymmetric 
relative to both axes of symmetry of the section. The maxi mum bending 
moment arising at the Juncture point of the contour for unit generalized 
deformation Y4 =1 is determined by the formula 

M= — (5.3) 

di | ( EJq ) + d 2 I ( EJ 2 ) 

in which and J 2 are the moments of inertia of the vertical and 

horizontal elements of the closed rectangular frame of width d-^, d 2 . 

If the shell consists of only rectangular plates having stiff Junctures 
5-, 3 S 2 3 

at the ribs J-, = -■*=—, J 0 = 

1 12 * ^ 12 


and therefore 



NACA TM 223k 


25 


M = 


E 

d^ 5-j^ + d 2 1 8g3 


(5.4) 


Expanding equations (2.1l) for the case of torsion of the shell, 
computing the coefficients of these equations "by formulas ( 2 . 13 )* and 
using the diagrams given in figures 13 * 14 , and 15 * we obtain 


gf Edl 2 d 2 2 (F 1 + F 2 X" - | G(d 2 2F 2 


+ d 2 2 F 1 )u 4 - | G(d 1 2 F 2 - d 2 2 F 1 > 1 * 


- | G(d 1 2 F 2 + dgXX* + P 4 = 0 

| G(d 1 2 F 2 - cL^FjlV + | <di% 2 + d 2 2 *i)V 

+ | G(d 1 % 2 - d 2 2 F 1 )Y 4 " + qi = 0 (5-5) 

| G(d 1 2 F 2 + d 2 2 F 3 )u 1 ,» + | G(di 2 F 2 - ch^iK" + | &(d 1 2 F 2 + d^X" 


96 

di | ( EJi ) + d2|(EJ 2 ) 


y 4 


04 = 0 


•where d^* dg are the widths of the vertical and horizontal plates, F-^, 
F 2 , the cross-section areas of these plates, P 4 = pj^z)* = c^(z )* 

04 = 04 ( 2 )* the generalized external forces referring to unit length of 
the shell and computed by the formulas 


P4 = / p(z,s)cp 4 (s) dSj 


o, = 


g(z,s>|f 1 (s) ds, % = J' o(z* s X(b) ds 


(5.6) 

in which p(z,s) and q(z,s) are the loads referred to unit area of the 
middle surface of the shell and acting in the direction of the generator 
and contour line respectively. 
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If the shell is acted upon hy a transverse vertical load q(z) 
applied in the plane of some vertical plate (fig. 12(a)) we have 


P X = 


0, 




(5.7) 


The equations ( 5 . 5 ) niay be represented in another form., hy t aking 
for the unknown function the longitudinal normal stress a = cr(z) and 
the transverse tending moment M = M(z) referring to the juncture 
point of the contour. We have 


V 






(5.8) 


Eliminating from equation (5-5) and with the aid of 

formulas (5-8) and denoting the angle of torsion of Yt (z) hy Q(z), 

qcU 

for p ] _ =0, qi = -<i k = — we obtain 


F 1 + F 2 „ 
a 


an 1 

*1*2 + 21-l 


8 

d l d 2 




- 8 


E 


^ 1*2 



E 1 
G- d g F q 


= 0 


(5.9) 


0" - 4 


(«i 


— V" V + — - — ^ = 

Gdg^P ± J 2GS' 1 d 2 


0 


The first two of these equations form a symmetric system relative to 
the functions a = cr(z) and M = M(z) and in the case of the homogeneous 
problem reduce to an equation of the fourth order for cr and M 

cr 17 - 2A 2 c" + B^cr = 0, M 17 - 2A 2 M" + B^M = 0 (5.10) 


or to an equation of the sixth order for the angle of torsion 6(z) 
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•where A and B are the generalized elastic characteristics determined 
by the formulas 


A 2 = 


2k 


|(EJi) + d 2 | (EJ 2 ) 


+ 

V»S*1 


Gdl 2 F 2 


(5.12) 


■Qk _ 384 1 

d^ 2 d2 (F-^ + Fp) dj + d^Jjg 


Equation (5.10) is the fundamental differential equation of the 
shell here considered having 1° of freedom for deplanation of the section 
and one for deformation of the contour. By determining from this 
equation the function ct(z) and then hy equations (5*9) the remaining 
functions M(z), 0(z), adding to the integrals of the homogeneous 

equations the particular solutions of the nonhomogeneous system ( 5 . 9 ) 
we obtain the general integral for the problem here considered. The 
general solution will be found with an accuracy up to 6 independent 
arbitrary constants. For determining the latter it is necessary in 
each particular case to use the boundary conditions which in the given 
problem will consist of three conditions on each of the ends of the 
shell. Thus for example, in the case of a shell hinged at the fixed 
ends z = 0, z = l the boundary conditions have the following form 

a = M = 9 = 0 for z = 0; 0 = M = 0 = 0 for z = l (5.13) 


Writing out these conditions we obtain six equations with six 
arbitrary constants of integration of the system (5 - 9 ) • 

We may remark that in addition to the accurate solution here 
described, equations (5»9) for the boundary conditions (5-13) can be 
integrated with the aid of trigonometric series. Bearing in mind 
conditions ( 5 * 13 ) the solution of equations ( 5 - 9 ) for q = constant 
may be written in the form 
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Z 

n=l,2. . , 


sin 


(2n - 1) 


JtZ 


M = sin 

n=l,2. . . 


(2n - 1 )itz 


(5.14) 


9 = / sin 


(2n - 1) 


JtZ 


n=l,2. . . 


n 


where o^, M^, and 9 n are the required expansion coefficients of 
equation (5.14). Substituting (5.14) in ( 5 . 9 ) we obtain 


. > g (a. - iftg _8_ Sa , 0 

6 l 2 Jtd^&i - 1) 


7T + B - 8 Hr- + -I 

d l d 2 b \ 1 1 \Gd 1 ^F 2 Gdg^F- 




(5.15) 


4e 

Gd^i 


Jt(2n - 1) 


= 0 


(2n - l)‘ 


n 


1 % 



gq 

GFid2Jt(2n — l) 


= 0 (n = 1,2,3,.. .) 
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From these equations and series (5»l4) the required functions c(z), 

M(z), and &(z) may he determined with any initi ally given degree of 

accuracy in the case where the shell is acted upon by a vertical 

uniformly distributed load q in the plane of any vertical plate 

(fig. 12). We may remark that the series (5.14) for equations (5*15) 

possess very good convergence and for practical purposes may be 

restricted to the two or three first terms of the expaasion, equation (5.14). 

Equations (5.9) represent a particular case of the general eight 
term equations proposed by the author as early as 1931 (references 2, 3, 4). 
These equations, thanks to the independent characteristics (elastic E, G-, 
geometric d^, dg, F-^, Fg, Jj, Jg ), permit considering practically 

important problems on the computation of thin-stalled aeronautical and 
structural shells of closed rectangular cross-section. Ulus, for example, 
if the shell is reinforced by ribs or diaphragms then for a sufficient 
stiffness of the transverse frames we can set EJ^ = EJg = equa- 
tions (5.9) then go over into the equations of a shell with rigid 
(nondef ormable ) contour. 

In the case of a shell consisting only of rectangular thin plates 
hinge-connected along the ribs, the second of equations (5*9), which 
represents according to the general theory given in the works of the 
author (references 3 and 4) the equation of continuity of the angular 
deformations on the contour line of the cross-section, drops out. In 
the remaining two equations the moments M must be set equal to zero. 

We obtain two equations for a and 0 referring to the shell with 
zero moment. 

If in equations (5*9) the elastic generalized characteristics GF^. 
and GFg referring to the deformations of shear are set equal to 
infinity, we shall have the case considered in detail in the work of 
the author (reference 3 ) and referring to a shell having the properties 
of deforming with longitudinal tension (compression) of the fibers of 
the shell and with bending of the elementary closed frame of the cross- 
section. From this short analysis the meaning of the generalized 

elastic characteristics A 2 , of one of the equations (5.10), (5.1l) 
is explained. It is clear that in the case of the loading of the shell 
by a nonsymmetrical load the stresses and deformations deter min ed on 
the basis of the solution here described must be combined corresponding 
to the stresses and deformations arising as a result of bending of the 
shell as a beam. The elementary theory of the bending of beams, as 
is seen from the above discussion, constitutes a particular case of the 
general theory here given. 


6. Example of the Computation of a Shell of Composite Section 

We consider a shell with the section shown in figure 18 . We assume 
that this shell is attached in some manner to the supports z = 0 , 
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z = l and is acted upon "by a single transverse load, perpendicular tc 
the axis of symmetry of the section. Of the required generalized 
displacements in this case of loading two will he longitudinal, Uq, 

Up, and three transverse, Yq, Yg, Y 3 . Figure 19 shows the anti symmetrical 

diagr ams of the orthogonal functions cpq, cp^ corresponding to the 
displacements Uq and Up. Of these functions q>q represents the law 

of linear distribution of the displacements over the section z = constant. 
The function represents the deplanation of the section. The 

generalized longitudinal forces corresponding to the generalized 
displacements Uq arid will he the tending moment about the axis 

of symmetry and the bimoment. 


Figure 20 shows the diagrams of the functions ijr^, and lfr^ 

corresponding to the generalized transverse displacements Y-^, Vp, 
and Y^- The magnitude Yq represents the deflection of the shell in 
the direction perpendicular to the axis of symmetry, Y2 — the angle 
rotation about the point c on the axis of symmetry, Y^ the generalized 
deformation of the contour of the section. The chosen functions 
i|fp, and ijf 3 are orthogonal. 


Expanding equations ( 2 . 1 l), computing the coefficients of these 
equations by formulas (2.13) and using the diagrams given for the 
functions cpq 3 <P 2 , ^PqS <£2** ^1* ^ 2 > ^3 assuming that the plates 

are hinge-connected to each other (the shell has zero moment) we obtain 
for the dimensions and loads shown in figure 18 and for E = 2 G the 
equations 

Fd%, " - 4 FU-, + H Up - 4 FY p * = 0 
6 x 25 c 


3 ^ Uq* + — Fd^Up" - FUp + — FYq* - Y 2 * + ^ FdV3" = 0 

25 1 75 2 62 2 25 1 8 2 * 3 


4 FU 1 * - — FUp* + 4 FYq" = -(Pq + Pp + P 3 ) i 


( 6 . 1 ) 


^ FdUp* + ~ Fd 2 V p " = 
o ^ 16 


3>/3 „ . 


5/3 


T dPi + ^dP 2 + ^dP 3 li 


^ FdU 2 ' + f- Fd 2 V 3 " 


|,p aljp.ii,')! 

,16 1 8 2 g 3 y g 
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•where F = d 6 , 5 — thickness of the plate, P^, P 2 , and the unit 

transverse loads. We take 8 = 0.01 m, d = 0.50 m and 2 = 3*00 m. 

We as sume that the shell is hinge -connected cn. the supports z = 0, 
z = 2. Then for z = 0 , z = 2 

Uj* = U 2 * = Y 1 = Y 2 = V 3 = 0 (6.2) 

For these boundary conditions equations (6.1 ) can be solved by the 
"approximate" method in the form of trigonometric series for TJ 
in cos (njtz/2) and for V in sin (nitz/2). Restricting ourselves to 
the first terms of the expansion, assuming P^ = constant, P 2 = constant, 

P 3 = constant and solving the system of algebraic equations which are 
obtained from equations ( 6 . 1 ) on substituting in them 


JJ ± = A x cos y-,..., V3 = B 3 sin P x = ^ sin (6.3) 

we shall have 

Ui = | (-616P-L - 213P 2 - 213P3) cos T" 

u 2 = 1 (l24Pq - 363P2 + I26P3) cos 3 S. 

Y ± = i (527Pq + 51^2 + 527P3) sin ^ ( 6.k ) 

V 2 = | (-318P-L + 162P 2 + 432P3) sin ^ 

Y 3 = | (163P X - 536P 2 + 245 P 3 ) Sin 
The normal stresses a(z,s) are determined by the formula 
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where cp^ = q>^(s) and q>> = cp 2 (s) are the generalized coordinates 

represented in figure 19- Substituting the obtained values U-^ = U-j^z) 
and U 2 = UgCz ) in equation ( 6 . 5 ) we obtain 

a 1 = (+^5Px + 223P 2 + 223P 3 )qpi(s) sin 

cr 2 = (-130P! + 380P 2 - 132P 3 )<p 2 (s) sin (6.6) 

From the first of these equations are determined the normal stresses o ^ 
referring to the bending of the shell as a beam (cpq represents the 
distance from the axis of symmetry of the section to the point of the 
contour.) The second of formulas ( 6 . 6 ) refers to the stresses arising 
as a result of the deplanaticn of the section. 

The maximum stresses Oj and c 2 arise in the middle transverse 
section ^for z = P-^ P 2 , P 3 are the intensities of the transverse 

loads in kilograms per meter. The stresses a determined by formulas (6.6) 
are obtained in kilograms per meter 2 . By comparing cr-p cr 2 it follows 

that the stresses due to the bimoment are of the same order aB the 

stresses cr-j_ due to the bending. Thus for example if the shell is 
acted upon by only a Bingle transverse load P 2 , then setting in 

equation ( 6 . 6 ) Pj = P 3 = 0 we obtain z = \ 

= 223P 2 <pq(s),, ct 2 = 380P 2 cp 2 (s) ( 6 . 7 ) 

For the lower point a of the intersection of the end vertical 
plate with the inclined one according to figure 19 we shall have 

c 1 a = +223P 2 d = 223 X 0.5P 2 = 1U.5P 2 

(in kg/m^) ( 6 . 8 ) 

a 2 a = - 38 op 2 ^ d = -380 x 0 .i4p 2 = -53-3P 2 

The total stress at this point will be 


cr 8 - = 58 . 2 P 2 (in kg/m.2) 
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It is seen from this example that the factor of deplanation of the 
section is of essential significance. We may remark that from, the 
formula (6.6) for erg it follows that there exists a system of loads P-^ 

P 2 , P.o for vhich deplanation of the section is absent. These loads in 
the' given example are evidently connected by the relation 


-130P], + 380P 2 - 132P 3 = 0 


(6.9) 


7* On the Principle of Saint— Tenant in the Theory of Shells 

Let the shell with the cross-section shown in figure 11 have the 
dimensions 


= 8 2 = 8 (7.1) 

in the absence of ribs assume the 
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where 6 is the thickness of the plates of the shell which in the 
given case we assume the same. Further let E = 2G. 

We consider an infinitely long shell and assume that at the initial 
section there is applied a system of longitudinal forces reducing to a 
bimoment and therefore representing the balanced longitudinal load corre- 
sponding to the deplanation of the section. The distribution of the 
stresses a due to such load statically equivalent to zero over the 
section z = constant is defined by the antisymmetrical diagram shown 
in figure 13 • We assume further that the system of longitudinal forces 
is such that the normal stress due to these forces at a juncture point 
of the cross— section is in absolute value equal to unity. The 
tangential stresses at the initial section z = constant will be considered 
equal to zero (fig. 21 ). 

Setting in equations (5*5) p = 2, p = q = 0, using the dimen- 
sions (7.1 )* integrating these equations and remembering the above 
boundary conditions at the intial section z = 0 and the condition of 
finiteness at z = °° we obtain a definite solution for the infinitely 
long shell of square section loaded at the initial section by a single 
balanced longitudinal load. 

In figure 22 are shown the graphs of the variation of the 
magnitudes -tj- proportional to the normal, stress a = a{z) as a 


(7.2) 


di — <±2 — dj 

The moments of inertia and 

values 
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function, of the position of the section along the shell. In figure 23 
are given the graphs of the variation of the magnitude Y^/d 2 propor- 
tional to the transverse bending moment M and characterizing the 
deformation of the contour of the shell. On these graphs along the 
axis of abscissas are laid off the relative coordinates £ = giving 

the distance in fractions of the width of the plate d from the origin 
z = 0 to the section z = constant under consideration. On the axis of 
ordinates are laid off in figure 22 the values u*/d and on figure 23 * 
Yj,/d 2 computed for different values of the ratio 8/d of the thickness 
of the plate to its width. 

From the graphs for a = Eu* shown in figure 22* it is seen that 
for a thickness of shell 8 = O.Old the normal stresses a due to 
the balanced longitudinal load do not have the character of local 
stresses and decrease very slowly with Increasing distance from the 
point of application of this load. Even at sections at a distance z = 8d 
from the initial section these stresses constituted 20 percent of the 
given stress at the section z = 0. The rate of decrease depends on 
the ratio 8/d. On increasing this ratio the rate of decrease increases. 
This increase however does not occur rapidly enough for the stresses 
dete rmin ed by the fourth term of formula ( 5 . 1 ) and, arising from the 
deplanation of the section to have the character of local stresses. 

For the case ^- = 0.1 the stresses a in the section z = 2.70d still 
d 

remain appreciable. 

From the example given it follows that the principle of Saint— 

Yenant, as already several times noted in our previous works (references 1 
and 4)* has a restricted field of application in the theory of shells. 
Thus* for example, for the shell here considered for ^ = 0.01 and for 

a length not exceeding 15 to 20 times the width of a single plate the 
stresses a due to the balanced longitudinal load applied at both 
ends of this shell, even at the middle cross-section at the greatest 
distance from the ends of the shell constitute 25 to 30 percent of the 
given stresses at the end sections z = 0, z = l. Thus, in practical 
problems the stresses due to the balanced longitudinal (or transverse) 
loan for thin shellB of finite length will not be local stresses. The 
principle of Saint— Yenant applicable in a practical sense to dense bodies 
can lead, on application to shells, as has been shown in our work on 
thin-walled rods, to false results. It follows that the elementary 
theory of the bending of shells as beams represented by the first three 
terms of the general law ( 5 . 1 ) and based essentially on the principle of 
Saint— Yenant iB applicable to particular cases, namely, when the external 
forces and the boundary conditions are such that the bimoment of the 
longitudinal forces associated with the deplanation of the section is 
equal to zero, that is in those cases where the system of equations (5*5) 
has a zero solution., In the light of our discussion above, the principle 
of Saint— Yenant has a restricted field of application also in the theory 
of shells. Thus if the closed circular cylindrical shell is subjected 
at the sections z = 0, z = l to normal stresses varying along the 
section according to the law cos^Q where 0 is the central angle, and 
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reducing therefore to a "balanced longitudinal load, the stresses in the 
center cross-section z = ^ for a length of shell 1 = 20E and thick- 
ness 5 = 0.01R (R is the radius of the circle) constitute, as shown 
"by the"'investigations of the author, from hO to 50 percent of the given 
stresses at the end sections. With decrease in the thickness 5 of the 
shell the stresses in the section z = ^ increase. In the case of a 


shell with zero moment the stresses over the entire length remain constant . 
In the cylindrical or prismatic shell the stresses due to the "balanced 
longitudinal forces at each of the sections z = 0, z = l have the 
character of local stresses only in the case of a sufficiently thick shell 

for > -zk , or if the thin shells for a thickness § < are reinforced 
K 3U 30 

over the length "by sufficiently rigid ribs. 


The problem, brought out above connected with the principle of 
Saint— Venant can be easily investigated with the aid of equations (5-5) 
by varying in these equations the magnitudes Jq, J 2 representing for 

a ribbed shell mean (referring to unit length) reduced moments of inertia 
of the longitudinal sections of the shell. With increase in the mean 
bending stiffnesses EJq and EJ^ of the elements of the transverse 

closed frame the rate of decrease of the stresses due to the longitudinal 
balanced load Increases. For EJq = EJ 2 = °°, that is in the case of a 

shell with rigid contour, the additional stresses connected with the 
deviation from the hypothesis of plane sections for longitudinal balanced 
load will have a local character. This remark refers to shells of 
closed profile. In the case of shells and thin-walled rods of open 
profile the stresses due to the bimoment, even in the absence of 
deformation of the contour of the section, extend over a considerable 
part of the length of the shell and do not have the character of "local" 
stresses . 


8. General Theory of the Vibrations of Prismatic 
Multiply-Connected Shells 

The general differential equations (2.1l) for the coefficients of 
these equations ( 2 . 13 ) and free terms ( 2 . 15 ) refer to the problem of the 
equilibrium of thin-walled spatial multiply— connected prismatic shell. 
These equations are generalized also to the theory of the vibration of 
such shells. For this it is necessary in equations (2.1l) to understand 
by the forces p^ and inertia forces arising from the vibrations 

of the shell for given forms, with an accuracy up to a certain number of 
parameters, of the longitudinal and transverse displacements. 

Considering the displacements u and v of any point of the 
surface of the shell as a function of three independent variables z, s, t 
where t is the time, 7 is the weight of unit volume of the material, 
and g the acceleration of gravity, we can write 
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P 


75 8^ 
S' st 2 ' 


76 8^1 

8 at 2 


( 8 . 1 ) 


Prom the above formulas are determined the surface Inertia forces, that 
is the forces referred to unit surface area. These forces act as 
follows: p" in the direction of the generator, and q in the direction 

of the tangent to the contour line of the tangent to the contour line 
of the cross-section. Setting, according to equations (l.l) and (l.2) 


u = U i (z,t)cp i (s) (i = 1,2,..., m), v = V k (z,t)iif k (s) (k = 1,2,..., m) 

( 8 . 2 ) 

Substituting these values in (8.l) and computing the free terms of the 
equations (2.1l) by formulas (2.15 ) we obtain 


( z,t ) = - ^ — ^U ± (z,t) J cpq(s) cpj(s) dF 

at 


q k (z,t) = - ^ ^2 J \( 8 )'M s ) ^ 


(8.3) 


Or taking into account formulas (2.13) 


(z,t) = — — / a.^ 1 U^(z,t) (i - 1,2,..., m) 

^ g at 2 


(8.4) 


q h (z,t)=-^> 1 ^. — V k (z,t) (k = 1,2,..., n) 

s at 

Equations (2.11 ) in the case of the dynamic problem will have the form 


I ZmiV - 21 - 1 Z - 2 


°jkV + 5 pj° - 0 


(8.5) 




c hi U 


i * + rh -t Vk " ~ S 2 a hk v k “ g r hk y k + g ^ = ° 
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■where the derivatives with respect to the coordinate z are denoted "by 
primes and those with respect to the time t by dots, pj° = Pj°(z,t), 

q k ° = q^ 0 ( z , t ) denote the generalized external longitudinal and 

transverse forces considered as functions of z and t and computed by 
formulas ( 2 . 15 ). If these forces are absent we have homogeneous equations 
in partial derivatives referring to the free vibrations of shells. These 
vibrations, as is seen from equations ( 8 . 5 ), will be harmonic. 

Using the method of the separation of variables we obtain a system 
of homogeneous ordinary linear symmetrical differential equations with a 
parameter determining the frequency of the natural vibrations of the 
shell. With these differential equations and the associated homogeneous 
statical, kinematic, or mixed boundary conditions (depending on the 
character of the problem), the frequency and the form of the vibrations 
of the multiply— connected shell are determined as of a spatial, 
discretely continuous system, the cross-sections of which do not remain 
plane. 

If the boundary conditions are such that for z = 0, z = l 

Uj* = 0 (i = 1 , 2 ,.../ m), Y k = 0 (k = 1 , 2 , . . ., n) ( 8 . 6 ) 

the fundamental functions of the homogeneous differential equations 
will be 


U lv (z,t) = Aq cos sin u^t, Y kv (z,t) = Bjj. sin sin o^t ( 8 . 7 ) 

Substituting ( 8 . 7 ) in equations ( 8 . 5 ) we obtain with Pj° = q k ° = 0 
for A^, Bjj. a symmetrical system of homogeneous differential algebraic 

equations. Equating the determinant of this system to zero we obtain the 
characteristic equation for the frequency of vibrations uiy. Since the 
initial equations are symmetrical all the frequencies of vibration will 
always have real values independently of the shape of contour of the 
shell cross-section. 

The investigations conducted by the author on the vibrations of 
thin-walled rods and shells show that the fundamental frequency of 
vibration will always be below that frequency which is obtained on the 
basis of the existing elementary theory of Rayleigh. In particular, 
for a shell of the type of a wing of an airplane, this frequency obtained 
on the basis of the theory here developed may be 3 to 4 times below the 
frequency obtained from the usual computation. 
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9. General Remarks 

The method here developed can be used to obtain the solution of the 
plane problem of the theory of elasticity for a rectangular region. Thus, 
for example, for the case of the plane parallel state of the rectangular 
plate the biharmonic problem leads to a system of symmetrically constructed 
ordinary differential equations 


21 V 


ffl. tp. dF 



+ 



vO f V<Pj ® 





dF - 





+ T lc"/ W* “ 

k 


(9.1) 



q h = ° 


(i,j = 1,2,..., m; k,h = 1,2,..., n) 

where v is the Poisson coefficient, = ^(z) (i - 1,2,... m), 

Yjj. = Y k (z ) (k = 1,2,... n) are the required generalized longitudinal 
and transverse displacements respectively, cp^ = cp^(s) (i = 1,2,..., M), 

^k = ^(s) (k = 1,2,..., n) are given functions depending on the 

position of the point along the width of the plates and determining 
according to equation (l.l) the state of the longitudinal and transverse 
displacements over the section z = constant. The functions cp^(s) 
and ikjj-C s ) 011 'the longitudinal edges of the plates must satisfy only 

the kinematic conditions. In choosing these functions the rectangular 
plate can be divided along its width by a number of narrow strips and 
within the limits of each such strip the assumption of linear distribution 
of the longitudinal and transverse displacements can be made within the 
limits of the width of a strip. In the particular case of m = 1, n = 1 
we obtain a system of two differential equations of the second order. 

These equations for = i)r - 8 and for the coordinate s measured 
from the center point of the width of the plate will refer to the case 
of the bending of a plate. If the coordinate s is measured from the 
lower edge of the plate, then for q>j_ = \|r = s we shall have equations 
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referring to a narrow rectangular plate which over its entire lower edge 
is rigidly fixed against vertical and horizontal displacements, etc. 

-We may remark, that, in contrast to equations (2.11 ), equations (9.l) 
were obtained with account taken of the tensile deformations of the plate 
over its width. Equations (9»l) in agreement with what was said above 
represent the generalized equations of the theory of bending of beams 
and make possible the solution of a number of problems. To such problems 
for example belong the contact problem on the computation of a beam 
lying cn an elastic base not following the hypothesis of Winkler, the 
problem of the local stresses in a narrow plate loaded by two concentrated 
forces applied on one vertical along the longitudinal edges of the plate 
and directed toward various sides, etc. 

For the cylindrical shell of arbitrary contour the differential 
equations in the variables u = u(z,s), v = v(z,s), in accordance with 
the assumptions here made for 5 = constant have the form 


i~i 5 2 u „ ci^u „ d 2 v 

E — — + G + G 


+ P = 0 


dz‘ 


ds c 


S zds 


■ (9.2) 


-.2 

<3z5s 


52 


0^I + 

c)z 2 12(1 - v 2 ) 





E §1 + <L_ I 

as a 8 e 


+ Q = 0 


where P and Q are terms depending on the surface load, E = E(s) is 
the radius of curvature. Equations (9.2) in contrast to the equations 
of Love have a symmetrical structure. It is not difficult to show that 
if the cylindrical shell is considered as prismatic and the number of 
sides is increased to infinity the system (2.1l) in the limit passes 
over into the equations (9*2). Equations (9-2) can be integrated by the 
method of the separation of variables. Setting (for P = Q = 0) 


u = Z(z)cp(s), v = Z*(z)iKs), Z" = X 2 z(z) (9.3) 


where X is an undetermined multiplier, we obtain 


EX^cp + Gcp" + Gt* = 0 


Gq>’ + GX 2 ^ + 


62 


'd d c 
vr E - 


1 d 
+ — 


12(1 - v 2 ) V dB “ as 2 T E de ^ 


d d t d ^ 

ds 2 ds ds E 


= 0 


(9.4) 


40 


NACA TM 1234 


The parameter X is found "by solving the correspond in g homogeneous 
“bo undar y problem described by the equations (9*4) and the homogeneous 
boundary conditions on the rectangular edges of the shell. In the case 
of a shell with rigid contour the problem reduces to the integration of 
the equations 


aCo&2 + j£rsl =Q 1 |£ _ I |L + S" x t + + 0"cd* = 0 

dz ob ^ os G dz 


(9-5) 


for the additional integral conditions 


'£0 

dz 2 


x dF + qx + 


dT 

dz 


= 0, F — — ^ y dF + q + 

J dz 2 7 


dT 1 

dJ y J " 


0 


(9.6) 


'£0 

6z 2 


a) dF + m + 


dT 

dz 


cu 


= 0 


In equations (9.5) | = £(z), q = t}(z), and 0 = 0(z) are respectively 

the deflections and angle of torsion, x = x(s), y = y(s), (0 = a>(s) are 
respectively the cartesian coordinates and the sectorial area. In 
equations (9.6) the integrals are taken over the entire contour of the 
section for the variable dF = Sds, T = T(z) denoting the shearing 
force. 


Translation by S. Reiss 
National Advisory Committee 
for Aeronautics 
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